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Abstract 

Entanglement witnesses are nonpositive Hermitian operators which can detect the 
presence of entanglement. In this paper, we provide a general parametrization for or- 
thonormal basis of C" and use it to construct projector-based witness operators for 
entanglement detection in the vicinity of pure bipartite states. Our method to pa- 
rameterize entanglement witnesses is operationally simple and could be used for doing 
symbolic and numerical calculations. As an example we use the method for detecting 
entanglement between an atom and the single mode of quantized field, described by the 
Jaynes-Cummings model. We also compare the detection of witnesses with the negativ- 
ity of the state, and show that in the vicinity of pure stats such constructed witnesses 
able to detect entanglement of the state. 
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1 Introduction 

The interest on quantum entanglement has dramatically increased over the last two decades 
due to the emerging field of quantum information theory. It turns out that quantum en- 
tanglement provides a fundamental potential resource for communication and information 
processing [1, 2, 3]. A pure quantum state of two or more subsystems is said to be entangled 
if it is not a product of states of each components. On the other hand, a bipartite mixed 
state p is said to be entangled if it can not be expressed as a convex combination of pure 
product states [4], otherwise, the state is separable or classically correlated. It is, therefore, 
of primary importance testing whether a given state is separable or entangled. For systems 
with dimensions 2 Cg) 2 or 2 (g) 3, there exists an operationally simple necessary and sufficient 
condition for separability, the so called Peres- Horodecki criterion [5, 6]. It indicates that a 
state p is separable if and only if the matrix obtained by partially transposing the density 
matrix p is still positive. However, in higher dimensional systems this is only a necessary 
condition; that is, there exist entangled states whose partial transpose is positive. 
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Peres-Horodecki criterion for separability leads to a natural computable measure of en- 
tanglement, called negativity [7, 8, 9]. Negativity is based on the trace norm of the partial 
transpose p Tl of the bipartite mixed state p, and measures the degree to which p Tl fails to be 
positive, i.e. the absolute value of the sum of negative eigenvalues of p Tl 

m s i£J-i, (1) 

where || p Tl || i denotes the trace norm of p Tl . Vidal and Werner [9] proved that the negativity 
N{p) is an entanglement monotone and therefore it is a good measure of entanglement. 

The most general approach for detecting entanglement is using entanglement witnesses 
[6, 10, 11, 12, 13]. Entanglement witnesses are operators that are designed directly for dis- 
tinguishing between separable and entangled states. By definition, we say that a Hermitian 
operator W defined on the product space T-L = T-L\ ® H2 is an entanglement witness if and 
only if: 1) Tr(Wa) > for all separable states a £ S, and 2) there exists at least one entan- 
gled state p such that Tr(Wp) < 0. The negative expectation value is hence a signature of 
entanglement, and for a state p with Tr(Wp) < we say that it is detected by W. It turns 
out that a state is entangled if and only if it is detected by some entanglement witnesses W 
[6]. 

An important class of entanglement witnesses is the so called projector-based witness. 
Given a pure entangled state \ip) £ Hi ® W2, its entanglement witness is given by 

^ = fc(V0i-|VXV>l, (2) 

where k(tp) comes from the maximal fidelity between and a product state, i.e. 

fc(V)=ma X |<e,/|V>>| 2 . (3) 
|e,/}es 

In this paper, we provide an explicit parametrization for the general orthonormal basis of 
the Hilbert space C™. Naturally, such a parametrization is closely related to the parametriza- 
tion of unitary matrices U{n) [14], and therefore such basis requires, in general, n 2 real 
parameters, i.e. the dimension of unitary group U{n). This parametrization can be useful in 
problems arising in quantum information theory. For instance, they can be used to construct 
maximally entangled states (or generalized Bell states) of a bipartite system, or Grecnberger- 
Hornc-Zeilingcr (GHZ) [15] states of a multiqubit system. We use these maximally entangled 
states and construct projector-based witnesses for detecting entanglement. Our method to 
construct entanglement witnesses is operationally simple, in the sense that they can be stored 
in a computer and that could be used for doing symbolic and numerical calculations. As 
an example we use the method for detecting entanglement between an atom and the single 
mode of quantized field, described by the Jaynes-Cummings model [16]. We also compare the 
detection of witnesses with the negativity of the state, and show that in the vicinity of pure 
stats such constructed witnesses able to detect entanglement of the state. 

The organization of the paper is as follows. In section 2 we provide an explicit parametriza- 
tion for the general orthonormal basis of the Hilbert space C". In section 3 we construct 
projector-based entanglement witnesses. In section 4 we start by reviewing the Jaynes- 
Cummings model and its solutions and calculate the negativity of the final state and use 
the constructed witnesses to check its separability. The paper is concluded in section 5 with 
a brief conclusion. 
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2 Parametrization of the general orthonormal basis 

The aim of this section is to introduce the most general orthonormal basis B for C™. We 
denote this basis by 

B = {|$ (0) >, |$ (2) ), • • • , l^™- 1 ))}. (4) 

Let {lei)}™^ 1 be the computational orthonormal basis for C". In this basis a general normal- 
ized vector |$(°)) can be expressed by 2n — 1 real parameters as (we do not remove the total 
phase) 



|$(°)) = cos^e^'leo) 

+ sin^cosflfe^V) 

+ ••• (5) 

+ sin^ • • • sintff 2 cos ft e^|e n _ 2 ) 

+ S m^...sinftsinfte i *»^|e n _ 1 ). 
Also we can find n — 1 vectors, which are orthonormal to the above state, as [14] 

l*Bi> - Jo) l*^ k«=^>=...=^ 1 =./2' fc = 1, ■ ■ ■ , » - 1, (6) 

where in the above formula one calculates first the derivative and afterwards the restriction 
to 7r/2. Despite the fact that the set {|<I>( )), |$fc+i)}fc=i ^ s orthonormal and constitutes a 
basis for C n , but they are not in general form, in the sense that any U(n — 1) transformation 
of the subset is also orthonormal to |<f>(°'). We therefore define the new vector 

l^ 1 )) as a linear combination of all vectors of the subset {|3>fc+i)}fc=i as below 



fe+i/ 

<&«) = cos^e^yW) 

+ sin^cos^e^Vi 05 ) 

+ ••• (7) 
+ sin^W • • • sin ft cos ft e^-3|<ft) 
+ sinfl^ • • • sin ft sin ft e**"-'!*^). 



Obviously, this vector is orthonormal to the |$(°)). Furthermore, there exist also n— 2 vectors 
orthonormal to l^ 1 )) as 

= Ji) l$(1)> ^^=...=e^ 1= „v fc = 1, ■ ■ ■ , » - 2. (8) 

By construction, these vectors are also orthonormal to the primary vector |$(°)). Again this 
new subset {|$fc+i)}fc=i is not unique and any U(n — 2) transformation of it, has also the 
same property. Therefore the third vector of the set B can be obtained by making the linear 
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combination of the above vectors as 

5(2) ^ 2> i*ai 



+ sin^ 2) cos^ 2) e^l*^) 



+ ••• (9) 
+ sin 9^ ••• sin ^ 4 cos 6™ 3 e*»-« 1 2 ) 
+ sin0< 2) • ■ -sin ft sin ft e^-3|<ft>. 



Taking the derivatives of |<1>( 2 )) with respect to 9^ for k = 1,- • • ,n — 3, and making linear 
combination of the obtained vectors, we get the vector |$( 3 '). Continuing this procedure, 
iteratively, we can find all elements of the orthonormal basis B, which in summary can be 
written as (m = 0, 1, • • • , n — 1) 



where 



+ sin^cos^e^l^™-^) 

+ ••• (10) 

+ sin<) • ■ -sineUcosei-i e^-^t^) 

+ sin^ • • • sineU sinei-i e^-'— |<™^i>. 

= 3^0 l<I>(m)> lfli=« a =-=e*-i=T/2. fc = 1, ■ ■ ■ , » - m - 1. (11) 

Here we have defined ^) = |efc_2). From this it is clear that for a given m, the number 
of parameters required to express |<£>( m )) as a linear combination of the vectors of the subset 

{1^1+1^)}!?=™ i s cc l ua l to 2(n — m) — 1, and consequently we need, in general, X)m=o(^( n — 
to) — 1) = n 2 parameters, i.e. n(n — l)/2 angles and n(n + l)/2 phases. This number is, 
actually, the dimension of the group of unitary transformation U(n). Indeed if we write such 
constructed basis as |<j>( m )) = ^™=o U m i\ei), (to = 0,1,--- ,n — 1), then it is this matrix 
U which is unitary with det U = Exp^z ^m=o Y^l=\ L _1 4™') • ^ follows therefore that the 

SU(n) transformation can be achieved if we add the requirement ^m^oS^Ti"" 1 = 0- 

For more illustration of the method we give below two simple examples. First let us 
consider n = 2. In this case we have 



|$ (0) ) = cos Pi e 4 * 



costfi e'^lei) + sin 0i e l01 |e 2 



= -sin6»i e (i</,0+?o) |ei) +COS6I! e^+^^ez). (12) 
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Also for n — 3 we have 

|$(°)) = cos6>i e'*°|ei) + sin0icos6> 2 e i01 |e 2 ) + sin 6>i sin 6> 2 c^ 2 |e 3 ), 

= -cos77isin6»ic t(0o+ « o) |ei) 

+ (cos r/i cos 6>i cos <9 2 e i(?0+ * l) - sin ??i sin 6> 2 e 4( « 1+0l) ) |e 2 ) 

+ (cos 771 cos 6»i sin 6> 2 e 4( « 0+<fc) + sin r/i cos 9 2 e l( « 1+<fc) ) |e 3 ), 

|$< 2 )) = siJUTisiJiflie^+^^lei) 

- (sin?7icos0icose' 2 e 4(Co+ « o+0l) + COS771 sin6» 2 c 4 (Co+«i+<M^) 

+ (-sin?7icos0isin6l 2 e l(Co+?0+ * 2) +cosr/iCOs6» 2 e 4(Co+ « 1+02) ) |e 3 ), (13) 

where for the sake of simplicity we have used 9^ = Oi, <f)f^ — (pi, 6^ = r)i, <j>^ = 4>^ = 
Co- 

3 Entanglement witnesses 

Now in this section we attempt to use the above parametrization to construct entanglement 
witnesses. Given a pure entangled state belongs to C" 1 ® C" 2 , its entanglement witness 
is given by 

^ = fc(V)l-|V)(V'|, (14) 
where k(tp) comes from the maximal fidelity between and a product state, i.e. 

fc(V) = ,ma X |(e,/|V)| 2 , (15) 

where S denotes the set of all separable states. This entanglement witness detects entangle- 
ment around the entangled state In general, it is not easy to calculate the constant k(ip), 
except for a two-qubit system that there exists a simple relation as 

where C(ip) is the so called Wootters concurrence [17], and for a pure state \ip) — a 00 |00) + 
a oi 1 01) + a 10 |10) + an |11), has the form C(ip) = 2|aoofln — aoi a iol- It is clear from equation 
(16) that k(ip) ranges from 1 to | as C(ip) goes from to 1, so that the minimum value for 
this constant happens whenever \i[>) is a maximally entangled state, i.e. a Bell state. On 
the other hand for a general bipartite system it is shown that the constant k(ip) equals to 
the square of the maximal Schmidt number of the state [18]. According to the Schmidt 
theorem, any bipartite pure state \tp) € C™ 1 <g> C" 2 can be written in the following form [19] 

n 

W = £v^k>h>i>, (17) 
i=i 

with 1 < n < min(ni,n 2 ) and Aj > with J27=i ^* = 1' an( I wnere an d {|i>i)}™ = i 

are the orthonormal eigenvectors of the reduced density operators p\ = Tr 2 (\4>)(ip\) and 
p 2 = Tri (\ip) , respectively. The number n is called the Schmidt rank of \ip). As a result 



Parametrization of projector-based witnesses 



6 



of the Schmidt rank, we can say that if a system has dimension m, then it can not be entangled 
with more than m orthogonal states of the another system. 

Motivated by this, we now use the orthonormal basis B and B' for the Hilbert spaces C ni 
and C™ 2 respectively, and define a bipartite pure state \^^) in Schmidt form as 

= cosai|$ (0) )|$' (0) > 

+ sina lC osa 2 |$ (1) )|$' (1) ) 

+ ••• (18) 
+ sinai • • • sina n _ 2 cos a n _i |$( n ~ 2 )) |$' ( ™~ 2) ) 
+ sinai • • • sina„_ 2 sina!„_i|<I> (n ~ 1) )|<l>' ( ™~ 1 ''). 
Obviously, the generalized Bell states obtained whenever all Schmidt numbers become 
which occurs when on — cos -1 . 1 . , , for i = 1, • • • , n — 1. We can, therefore, define the 

VII- 4+1 

witness operator based on the pure state |\I>M} as 

W 9W = fc(* w )l - |*W)^N| ; (19) 

where k(^ n ^) is equal to the square of the maximal Schmidt number of the ket |\E»[™]} [18]. 
Equation (15) guarantees that Tr (W^ w cr) = fc(*N) _ (^N| CT |^W) > o for all separable 
states (7 e «S. On the other hand since the Schmidt rank of every bipartite entangled pure 
state is greater than one, so all Schmidt numbers of the entangled bipartite pure states are 
less than 1 and therefore Tr (W $l „i |*N)(*W|) = fc(*N) - 1 < for every entangled state 
I'J'I™]), i.e. W^h detects entanglement of the pure state | 

Now let us consider the state p acting on the Hilbert space C™ 1 ® C™ 2 . Equation (19) guar- 
antees that if (Wy ln] ) p = Tr(W 4WP ) = fc(*N) - F 9W (p) < 0, where F vW (p) - 
is the fidelity between two states \^^) and p, then p is not separable and has some entan- 
glement. Therefore for a given state p acting on the Hilbert space C" 1 (g> C" 2 , we say that 
Wq,[ n ] detects entanglement of p if and only if the fidelity between and p is greater than 

fc($N), otherwise p is unentangled or its entanglement can not be detected by Wq,[ n \. 



4 Witnessing entanglement of the Jaynes-Cummings 
model 

In this section we use such constructed witnesses to detect entanglement of the Jaynes- 
Cummings model. The Jaynes-Cummings Hamiltonian between a two-level atom A and a 
single-mode quantized radiation field F is described by 

H = ^hojA^z + fc)ffl'o + %(cr+ <8> a + a- ® af). (20) 

This Hamiltonian acts on the product Hilbert space H A ®'H F . Here g is the atom-field 
coupling constant, uja = {e e — € g)/h is the atomic transition frequency, and ujp denotes the 
field frequency. The atomic "spin-flip" operators a + = \e)(g\, cr_ = |.g)(e|, and the atomic 
inversion operator a z = \e)(e\ — \g)(g\ act on the atom Hilbert space 1~L A — C 2 spanned by 
the excited state |e) — > (1,0) T and the ground state \g) — > (0, 1) T . The field annihilation and 
creation operators a and a) satisfy the commutation relation [a, a)} = 1 and act on the field 

Hilbert space V. F spanned by the photon-number states {|n) = |0)}^ . In the rest of 
this section, we consider the solutions of this Hamiltonian in two cases. 
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4-1 Case 1 

We first assume that the atom is initially prepared in the excited state p (0) = |e)(e|, and 
the field is initially in the number state p F (0) = \n)(n\. We also consider the effect of pure 
phase dccoherence on the Jaynes-Cummings model. In this situation, the master equation 
governing the time evolution for the system under the Markovian approximation is given by 
[20] 

ig = -i[H,p]-l[H,[H,p]], (21) 

where 7 is the phase decoherence coefficient. The formal solution of this equation can be 
expressed as 

p(t) = J2^-M k (t)p(0)M^ k (t), (22) 

fc=0 

where p(0) = p A (0) ® p F (0) is the initial state of the system and M k {t) is defined by 

M k (t) = H k cxp(-iHt) cxp(^iJ 2 ). (23) 
Then the time evolution of the system reads [21] 

p n (t) = £ n |e)(e|<g>|n)(n|+G„|e)(s|<g>|n)(n + l| 

+ G* n \g){e\® \n+l){n\+F n \g){g\® |n+l)(n+l|, (24) 

where we have defined 

A 2 

+ (2 - — — ) cos2Q„i exp(-2 7 ^ 2 
(2-2 cos 2Q, n t cxp(-2 7 tfi 2 )) , 



E n 


1 ( A 2 

" n 2+ 2o 2 


F n 


l<? 2 (n+l) 
4 tt 2 


G n 


aV( n + 1) | 

4ft„ 1 


where A 





g ^lo + ~ f(7~ - cos2n nt cxp(-27^ 2 )) + 2isin20„i cxp(-2 7 ift 2 ; 



I A 2 

as 



U{Pn) = g ^m fT \ (1 ~ cx P(- 2 ^n) cos 2n " f ) 2 + 4 exp(-4 7 in 2 ) sin 2 2Q n t ) . (25) 



Now we attempt to characterize the separability of the above state with witnesses defined in 
equation (19). To this aim we first consider a Bell state as 

|*P1) = _L (|$(o))|$'(°>) + |$(D)|$'«)) , (26) 

where {|<I>( )), l^*- 1 ^)} are two orthonormal vectors in the field space spanned by {|n), |n+l)}, 
and {|<&'^), |$' < - 1 ^)} arc two orthonormal vectors in the atomic space spanned by {|e), \g)}. 
The parametrization of both sets of vectors are given by equation (12). We find for the fidelity 
between p n (t) and |*[ 2 1) 

F 9 m (p n (t)) = \ + \Y 2 \ 2 F n ) + Rc (y*Y 2 S^-^+^-^q^ ; (27) 
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Figure 1: Negativity Af(px(t)) (upper panel) and the fidelity F^p] (pi(i)) (lower panel) are 
plotted as a function of t with 5 = 1,7 = 0.3, A = 1. The horizontal line in the lower figure 
shows the minimum fidelity that above it the state is detected to be entangled. 



where Y\ and Y% are defined by 

Yi = cos 6»i cos 9[ + sin 6 1 sin 9[ e^ 0+i '°\ 
Y 2 = sin 0i sin 0[ + cos 6 X cos 6[ e ^°+*o) . 

Now in order to get maximal fidelity, we can use the software MATHEMATICA, and maximize 
F^m (Pn(t)) with respect to all parameters of the state l^ 2 ]). In figure (1) we have plotted 
negativity J\f(pi(t)) (upper panel) and the maximum fidelity F 9 [a](pi(t)) (lower panel) as a 
function of t with g = 1, 7 = 0.3, A = 1. It is clear from the figure that whenever negativity 
is nonzero, the fidelity F 9 [2](pi(t)) is also greater than 1/2, which indicates the detection of 
entanglement by W^pj . 



4-2 Case 2 

Since, in practice, it is difficult to realize an atom in a pure state, therefore we now suppose 
that the atom is prepared, initially, in a general mixed state with the diagonal representation 

p A (0) = X\g){g\ + (l-X)\e){e\, AG [0,1], (28) 



but, however, the field is in a pure number state 

p F (0) = \n)(n\. 



(29) 
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Figure 2: Negativity Af(pi(t)) is plotted as a function of t with g = 1, A = 1. . 



However, we ignore the effect of decoherence and assume that the evolution of the system is 
unitary. Accordingly, the final state of the system can be obtained as [22, 23] 

p n (t) = A(K| 2 |e)(e|® \n-l){n-l\+A n B* n \e){g\® |n-l)(n| 

+ A* n B n \g)(e\ ® |n)(n- 1| + |B„| 2 | 5 >(g| ® |n)(n|) 

+ (1 - A) (|C n | 2 |e)(e| ® |n)<»| + C„£>*|e}< 5 | ® |n)(n + 1| 

+ C;£>„| ff )(e|(g) |n + l)(n| + |A>| 2 | 5 )(s|® |n +l)(n+ 1|) , (30) 

where we have defined 

C n = e- iUF ( n+1 / 2 )* ('cosfint-i^-sinOntV 

D n = _ ie — ("+ 1 /2)*^±I sin o„t, (31) 

where the Rabi frequency f2„, and the detuning parameter A are the same as before. For the 
above state the negativity can be expressed as 



Af{ P n{t)) = ~ (aA 2 |B„| 4 + 4(1 - A) 2 |C„| 2 |£>„| 2 - A|B„| S 

+ \ (v^l - A) 2 |C„| 4 + 4A 2 |A l | 2 |£? n | 2 _ (1 _ A)|C„| 2 ) . (32) 

In order to show the effect of mixing parameter A on the negativity of the system, we plot 
the negativity N(p n {t)) as a function of t and A with g = 1, A = 5 in figure (2). It is clear 
from this figure that when purity of this state is decreased, the negativity is also decreased. 
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t 



Figure 3: Negativity Af(pi(t)) (upper panel) and the fidelity F^in (pi(t)) (lower panel) are 
plotted as a function of t with g = 1, A = 5 and A = 0. 

Now in order to characterize entanglement of this state by using the notion of entanglement 
witness, we define l^ 2 ') as equation (26), but here {|$(°)), |$^^)} are two orthonormal vectors 
in the field space spanned by { | n~ 1), |n), \n + 1)} and {|$' ( ' ' ) ), l^'' 1 ')} are two orthonormal 
vectors in the atomic space spanned by {|e), \g)}. The parametrization of these two sets 
of vectors are given by equations (13) and (12), respectively. We obtain for the fidelity the 
following relation 

F^(p n (t)) = AQ(|X 1 | 2 |^|V|X 2 | 2 |B„| 2 )+Re(x*X 2 e^ 1 -^-^)A n 5;)^ (33) 
+ (1 - A) Q (\X 3 \ 2 \C n \ 2 + \X 4 \ 2 \D n \ 2 ) +Re(x*X 4 e 4 ^-^+^-^)C nJ D:)) > 

where 

Xi = cos 6»i cos 6' x + cos rj sin 6\ sin Q'^ ^ , 

X2 = sin 9i cos 6*2 sin 9[ + (cos 77 cos Q\ cos 2 e^° — sin 77 sin ^e 1 ^ 1 ) cos ^e 1 ^ , 
A3 = sin 6\ cos 82 cos & x — (cos ij cos #1 cos 2 e^° — sin ij sin ^e 1 ^ 1 ) sin 6' 1 e t ^ , 
A 4 = sin 0i sin 6 2 sin ^ + (cos 77 cos #1 sin # 2 e^° + sin r\ cos 2 e lfl ) cos Q'^' . (34) 

Again, by using the MATHEMATICA, we can maximize -F^pj (p n (t)) with respect to all 
parameters of the state IvE^ 2 ]). In figure (3) we have plotted the negativity Af(pi(t)) (upper 
panel) and the fidelity F^pj {pi(t)) (lower panel) as a function of t with g = 1, A = 5 and 
A = 0. In this case the state of the system is pure and it is clear from the figure that W^m 
enable to detect the entanglement of the system. On the other hand by increasing A from 
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Figure 4: Negativity Af(pi(t)) (upper panel) and the fidelity F^p] (pi(i)) (lower panel) are 
plotted as a function of t with g = 1, A = 5 and A = 0.2. 



to 1/2, the purity of the state is decreased, and the ability of W^pj to detect entanglement is 
decreased. For instance, we plot in figure (4) the negativity and the fidelity as a function of t 
for parameters the same as figure (3) but A = 0.2. We see that in this case there exist some 
entangled states that can not be detected by W^p]. 



5 Conclusion 

We have presented a general parametrization for orthonormal basis of C™. This parametriza- 
tion can be used to construct projector-based witness operators for entanglement detection 
in the vicinity of pure multipartite states. As an example we have used the method for de- 
tecting entanglement between an atom and the single mode of quantized field, described by 
the Jaynes-Cummings model. We have also compared the detection of witnesses with the 
negativity of the state, and have shown that in the vicinity of pure stats such constructed 
witnesses able to detect entanglement of the state. 
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